A new non-perturbative approach to Quantum Brownian Motion 
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Starting from the Caldeira-Leggett (CL) model we derive the equation describing the Quantum 
Brownian motion, which has been originally proposed by Dekker 0] purely from phenomenological 
basis containing extra anomalous diffusion terms. Explicit analytical expressions for the temperature 
dependence of the diffusion constants are derived. At high temperatures, additional momentum 
diffusion terms are suppressed and classical Langivin equation can be recovered and at the same 
time positivity of the density matrix(DM) is satisfied. At low temperatures, the diffusion constants 
have a finite positive value, however, below a certain critical temperature, the Master Equation(ME) 
does not satisfy the positivity condition as proposed by Dekker. 



The problem of Quantum Brownian Motion is a long- 
standing and challenging problem 0, [1] j as it forms the 
underlying basis of non-equilibrium phenomena such as 
dissipative and relaxation dynamics of quantum systems 
0i @|- Quantum dissipative dynamics has application in 
a wide variety of problems, starting from quantum cos- 
mological models [7, 8] to reaction rate theory Unlike 
the classical nonequilibrium problem, quantum dynamics 
has additional complexity due to the Heisenberg Uncer- 
tainty Principle. Numerous attempts have been made 
to address this problem, in a variety of ways. These 
include semiclassical approaches using Wigner distribu- 
tion function IQl , phenomenological models [1, [HI, , 
and using Boltzman's collision terms [l^]- On the other 
hand, there have been attempts to obtain the Master 
Equation (ME) describing the time evolution of DM of 
open quantum systems through toy microscopic models 
Il4| . where the environment is modelled by nonin- 
teracting particles or fields in thermal equilibrium. The 
ME is obtained for the reduced density matrix(DM) of 
the system by integrating out the environment degrees 
of freedom. New form of ME's are obtained from micro- 
scopic models where the diffusion constants have explicit 
time dependence 0, El]- Independently, a class of ME's 
have been put forward, [3] purely from mathematical 
consideration, which guarantee the positivity of the den- 
sity matrix. In other words, the time evolution of the 
system takes place only through the physical states, and 
the ME is commonly known to belong to the Lindblad 
class. Most of the MEs describing the non-equilibrium 
quantum dynamics suffers from the fact that either they 
do not belong to Lindblad class or correct classical limit 
can not be recovered at high temperatures. The ME 
derived in the original paper by CLpj], does not satisfy 
the positivity condition at high temperatures, although 
it recovers the classical Langevin dynamics at high tem- 
peratures. 

Independently a phenomenological model of quantum 
Brownian dynamics has been proposed by Dekker 
which contains extra terms describing momentum dif- 
fusion apart from ordinary diffusion, and the diffusion 
constants have to satisfy some conditions in order to pre- 



serve the Heisenberg's uncertainty condition as well as 
positivity of the DM 17, [l8|. However Dekker's model 



of quantum Brownian motion is purely phenomenological 
and the microscopic origin of the extra momentum dif- 
fusion constants are not known. Moreover, most of the 
derivations of the ME's rely on the fact that the auto- 
correlation function of the random force on the Brownian 
particle is short ranged at high temperature and that the 
system is assumed to be Markovian. But in practice, the 
effective action obtained from CL model has a memory 
kernel which decays as a power law at low temperatures 
[l^ . In fact, the random force in Quantum Brownian 
Motion is an operator and the force-force autocorrela- 
tion function has been evaluated exactly, which shows 
clear deviation from the Markovian limit '2(^ . This non- 
locality makes the problem very complex. In this letter 
we take a new approach, to analytically derive Dekker's 
form of quantum Brownian motion, starting from the CL 
model, where the analytical expressions of the diffusion 
constants are obtained for all temperatures. We show 
that at high temperatures, the ME satisfies the positivity 
condition and at the same time reproduces the classical 
Langevin dynamics. Although, the diffusion constants 
are finite and well behaved at all temperatures, the pos- 
itivity of DM breaks down below a certain temperature 
depending on the damping rate of the system. This is 
not very unexpected, since we know that at low tempera- 
tures, the long range memory effect becomes very impor- 
tant. This form of quantum Brownian motion is applica- 
ble to describe the relaxation phenomena and nonequilib- 
rium evolution of quantum systems, even at sufficiently 
low temperatures. We propose a new method to tackle 
the problem of nonlocality, using the equation of motion 
of the canonical coordinates. In this letter, we mainly 
focus on the quantum brownian dynamics in free space, 
however, our technique can, in principle, be extended to 
confined quantum particles. 

A very well known description for the dissipative phe- 
nomena and relaxation dynamics of classical system is 
given by the Langevin equation : 



q{t) = -jq{t) + eit) 



(1) 
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where, q{t) is the position of the Brownian particle, the 
dot denotes derivative with respect to time, 7 is the 
damping constant and 9{t) is related to the fluctuating 
force acting on the particle, whose autocorrelation func- 
tion is given by. 



{e{t)9{t')) =TS{t~t') 



(2) 



where F and 7 are related by the Fluctuation-Dissipation 
Theorem, which states that F = T M7, where T is 
the temperature of the bath and M is the mass of the 
Brownian particle. However, in most cases, the micro- 
scopic details behind this dissipation are not well known. 

The simplest microscopic model which describes the 
dissipative motion of a heavy particle in presence of a 
heat bath was first put forward by Ford, Kac and Mazur 
211 and later on applied to quantum systems by Caldeira 
and Legget [ij. In this model, full Hamiltonian has three 
different parts, system Hamiltonian Ha, the Hamiltonian 
of the heat bath Hb and the Hamiltonian Hj describing 
the linear coupling between the system and the bath, as 
given below 



Ha 
Hi 



E 



p? 
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(3) 
(4) 
(5) 



where, q and p are the position and momentum oper- 
ators of the heavy particle having mass M, V{q) is the 
potential in which the heavy particle is moving, Qi and 
Pi are the position and momentum operators of the bath 
oscillators whose mass and frequency of oscillation are 
given by and uji respectively and Ci is the coupling 
strength between the heavy mass and the z*'* oscillator 
in the bath. Following the work of Ford et. al., the time 
evolution of the position operator can be written as. 



Mq+ / dt'a{t-t')q{t') + V'{q)^ F{t) (6) 



where the prime denotes the derivative with respect to 
the position variable. The operator valued random force, 
F{t) is related to the statistical distribution of the bath 
variables Qi and PiS. In order to recover the classical 
Brownian dynamics with short time memory effect at 
high temperatures, the frequency dependence of the cou- 
pling is typically taken to be Cf (lu) / {niiUjf) oc M7, where 
7 is the damping constant, and the summation is then 
replaced by an integral over uj. According the Ehrenfest 
theorem, if the position and momentum operators are 
replaced by their classical equivalents, then Eq. [6] trans- 
lates to the classical Langevin Equation in the limit of h 
—f 0, where the force autocorrelation function takes the 



form of Eq. [2] with F = 2M-fkBT^ obeying classical 
fluctuation-dissipation theorem. 

The essential difference between the classical Langevin 
equation and the Quantum Master Equation is the fact 
that in the classical case, the fluctuations in canonical 
coordinates are controlled only by the scales set by the 
temperature and the dissipation constant. Whereas, in 
the Quantum case, there is an additional scale given by h 
which appears due the Heisenberg Uncertainty Principle. 
Moreover, the evolution of the operators themselves have 
to be unitary. Apart from the semiclassical techniques, 
new methods have been proposed to take into account 
the constraints imposed by 'uncertainty principle' in the 
quantum casejS^I. Among these methods, ME of the 
density matrix p is more suitable to describe dissipative 
dynamics of quantum systems. The quantum dynam- 
ics which is governed by the time dependent Schrodinger 
equation describes the pure state, whereas the dissipa- 
tive dynamics introduces the concept mixed state, where 
trace of is less than unity. The evolution equation of 
the reduced density matrix of the system can be obtained 
by tracing out the bath degrees of freedom. 

Using the Feynman- Vernon method [23| . the dissipa- 
tive dynamics of the reduced DM can be written in 
terms of the influence functional, which can be obtained 
after integrating out the bath degrees of freedom Qi. 
The main assumption of this method is that the sub- 
system is uncorrelated with the bath at the initial time. 
Hence the total density matrix at time t = is given by 
Pt(0) — Pa(0)®P_b(0), where pA, Pb and pt are the den- 
sity matrix of the subsystem, bath and the total system 
respectively. However, several authors have considered 
a correlated initial state, where the subsystem and bath 
density matrix cannot be factorized. For example. Hakim 
and Ambegaokar have compared the two cases of 
uncorrelated and correlated initial conditions and shown 
that, for the correlated initial condition, different tran- 
sient behaviours can be obtained at time scales larger 
than the inverse cutoff frequency of the bath, in contrast 
to the uncorrelated initial conditions. However, here we 
are following the usual procedure of Feynman and Vernon 
and we will show finally that we can obtain a consistent 
ME, describing the time evolution of the DM. The time 
evolution of the reduced DM is given by, 

P{qi,q2,t) ^ J dq[dq2J{qi,q2,t;q[,q'^,0)p{q[,q'^,0) (7) 

where the propagator J(gi, q2, t; q[, q'2, 0) is given by, 

J{qi,q2,t;q[,q'2,0) = J J VqiVq2Cxp Qs'off [gi, 92] 

(8) 

After integrating out the bath degrees of freedom, one 
obtains the nonlocal effective action corresponding to the 
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dissipative system: 



^^eff = ^iSA[qi] - SA[q2]) + j\^R + ^l)dT (9) 



where Sa is the action corresponding to the non- 
interacting system, and, 

^ -j^ j [q-{T)aR{T - s)q^{s)]ds 

^ -I I [q-{T)ai{T ~ s)q+{s)]ds (10) 



where, g± = gi ± q2 and the memory kernels are, 

otRiT) = coth ( I cosfo^ir) (11) 

a/W = -y]:5-^sin(w,T). 

i 

It is assumed that the oscillator frequencies are con- 
tinuously distributed from zero to a maximum frequency 
Uc, value of which depends on specific physical sys- 
tem. Cutoff frequency of the heat bath is chosen in 
such way that the characteristic time scale of the dy- 
namics of heavy particle is much larger than the coUi- 
sional time scale l/wc- Thus, the summations in Eq. 
1121 can be replaced by integrals by introducing an ap- 
propriate density of state, which will make the mem- 
ory function anaytically tractible. Hence, can be re- 
placed y, / dujF{uj), where F{uj) is the density of states 
of the bath. We choose a smooth Drude form for the 
density of states given by F{uj) — a;^/(ti>^ + uj-^) and 
F{Lo)C\Lo)/{2mio^) = (2M7/7r)(w2/(a;2 ^cj2)). At any 
non-zero temperature, using the Drude form of density of 
states the memory kernels can be evaluated analytically 
by using contour integrals in complex u plane and are 
given by: 



aRir) = M7w^ [cot(x) exp(-WcT) 



rnr/x 



-y 

^ ^1 (Wx)^ 



1 



exp [-{mr/x)(^cT] 



ai{T) 



MjujcTT- exp[-WcT] 

OT 



(12) 



(13) 



where x — hLUc/^ksT. From the above expression of an 
we can clearly see the emergence of two time scales - 
the first being the microscopic collision time, given by 
1/lOc and the other time scale given by temperature i.e. 
h/ksT. It is interesting to note that in the high temper- 
ature regime, an is short ranged over the thermal time 
scale while the timescale over which aj decays is given 
only by the collision time scale, which is independent of 
temperature. However, as the temperature is reduced, 
at some point, the thermal time scale dominates over 



the collision time scale, giving rise to a nonlocal memory 
kernel which has a power law decay. 

We now proceed to simplify the nonlocal action, taking 
advantage of the above mentioned short-ranged kernel. 
The dynamics over a timescale larger than ft/ZcsT can 
be studied by assuming smooth classical trajectories and 
expanding the dynamical variables in a Taylor series. 



9±(s) 



/=0 



II 



(14) 



where, g'-'-' denotes the Z*'* derivative of q with respect 
to time. We now insert Eq. [T^] and Eq. [T31 into Eq. 
[TOl Neglecting the total derivative terms which generate 
the boundary terms in the action, thereby leaving the 
Lagrangian invariant, we obtain, 



1=0 



i^'aR{T)df (15) 



The transient term in the integral comes in the form of 
exp{—{kBT / fijT . Thus, for r ^ h/ksT, we obtain. 



Sfl(r) = X 

fix 



eM-2x:(,^>)><^x:g)""<(2,„) 



-.(n„(-l) 

= 1 m=0 

Similarly, one can evaluate S/ which gives. 



S7(r) 



-q-{T)q+{T) 



(16) 



(17) 



where, the higher order terms in the expansion of q+{s) 
have been neglected, since they fall off as l/ujc- 

We have thus transformed a highly nonlocal action into 
a local action containing higher derivative terms. While 
in general, this may seem to be a very complicated ex- 
pression to work with, we show below that in at least the 
case of a particle in a harmonic potential attached to a 
heat bath and a Brownian particle there is considerable 
simplification of the expression which allows us to obtain 
exact results. 

For the free particle, general dynamical equation for 
any order of time derivative of the position coordinate 
can be written as, 



M ^ 

Substituting Eq. [18] in Eq. [TH we obtain, 



(18) 



Sfl(r) 
M-foJc 

hx 



MjUc 2 



hx 

q-\r) 



qt{r) 



(g)coth(g) 



(19) 
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Thus, the effective action is given by, 



n 



—q+q.~ -fMq^q+ 



ft2 



dT[q^_ 



aq_ 



(20) 



where a = *^°*^(2r^) ~ l]/7 ' assuming 7 < Wc- 

It is important to note that by using the dynamical equa- 
tion of canonical coordinates and resummation method 
we have converted the non-local action into an effective 
action which is local and quadratic in canonical coordi- 
nates as well as independent of the cut-off frequency coc 
of bath oscillators. Inserting the above form of the effec- 
tive action into Eq. [5] the time evolution of the DM can 
be evaluated. 

Following the prescription of CL we consider the 
change of density matrix from t to t + e within a small 
time interval e, in order to obtain the ME in differential 
form. To do so we expand both side of Eq. [7]upto leading 
order in e. Within the small time e, we approximate 
qi^ {x- x')/e = /3i/e and ^2 = (y - y')/e = /32/e- Now 
Eq. [7] reads: 



dj3+d(3-exp 



—M(3+{(3^ - 2-fX-e) 
he 



— {x_ -f a(}_) 



p(x-/3i,j/-/32,t)(21) 



where /3+ — (/3i + /92)/2 and /3_ = Pi - ^2- We expand 
p upto second order in /3, which is equivalent to expand- 
ing upto first order in e. After doing some algebra and 
performing the gaussian integrals we obtain the ME de- 
scribing the time evolution of the density matrix: 



dp 
di 



ih 
2M 
2i 



9^ 

dx^ 



9V 

,dp 



' dx 



dp 



dx dy 
d d 



+ -T-Dpq{x^y){^ + ^)+D,g{— + —Yp 



dx dy 



yfp. 



where, different diffusion constants are given by: 



Dpq — 



2kBT 
M7 

4fcBT 



( J:^) coth ( 

\2kBTj \2kBT 



- 1 



2kBT 



coth 



2fcsT 



Dpp = 2MjkBT 



, ^7 
'knT 



coth 



2kBT 



-3 



(22) 

(23) 
(24) 
(25) 



The above form of the ME has been proposed earlier phe- 
nomenologically 0] and the diffusion constants have to 
satisfy some conditions, in order to preserve the Heisen- 
berg Uncertainty Principle in quantum dissipative sys- 
tems 17|. This form of the ME can also be recast into 



the Lindblad form [l8| which maintains the positivity of 
DM during time evolution. Commonly known positivity 
condition is described by[l3, : 



A — DppD qq 



pq 



> 



(26) 



At high temperatures, the diffusion constants behave as, 
1 



Dpp = 2M-fkBT 



2 ( Jn_ 

3 \kBT 



Dqq = 



and 



Dnn ~ rp . The forms of these diffusion constants are 
very similar to those obtained earlier However, we 
would like to point out that the diffusion constant Dpq, 
in our case, is independent of the cutoff parameter Wc, 
unlike the earlier result. It is also interesting to note 
that the diffusion constants Dqq and Dpp both originate 
from quantum effects and vanish at high temperatures as 
l/kBT. However, at high temperatures, A approaches a 
value h?^'^ /12 which preserves the positivity condition of 
the ME. At zero temperature, the diffusion constants are 
finite, positive and proportional to h, similar to the form 
proposed by Dekker [2]. However, A becomes negative at 
zero temperature, violating the positivity criterion. We 
estimate a critical temperature Tq ^ 0.2h'y, below which 
above form of ME is not valid and transient behaviour 
as well as long range memory effects become crucial. 

To summarise, within our scheme of calculation, we 
have shown that the ME obtained from the microscopic 
CL model satisfies the positivity condition and belong to 
the Lindblad class above a certain temperature. Analyt- 
ical form for the diffusion constants have been obtained 
for any arbitrary temperature and are independent of the 
cutoff frequency of the heat hath. It is interesting to note 
that at high temperatures all the anomalous diffusion 
constants vanish as 1/T, which preserves the structure 
of classical Brownian motion. At the same time, the dif- 
fusion constants conspire in such a way that they satisfy 
the Dekker criterion A > 0, hence the time evolution of 
the open quantum system takes place only through phys- 
ical states. At low temperature, the diffusion constants 
have their origin from purely quantum effects. However, 
below the temperature ~ fi.7, Dekker's positivity condi- 
tion is violated, indicating that the long range memory 
effect may become important in the time evolution of the 
system. 
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